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We consider the quantum mechanical propagator for a particle moving in a d-dimensional Lorentz 
gas, with fixed, hard sphere scatterers. To evaluate this propagator in the semi-classical region, and 
for times less than the Ehrenfest time, we express its effect on an initial Gaussian wave packet in 
terms of quantities analogous to those used to describe the exponential separation of trajectories in 
the classical version of this system. This result relates the spread of the wave packet to the rate 
of separation of classical trajectories, characterized by positive Lyapunov exponents. We consider 
applications of these results, first to illustrate the behavior of the wave-packet auto-correlation 
functions for wave packets on periodic orbits. The auto-correlation function can be related to the 
fidelity, or Loschmidt echo, for the special case that the perturbation is a small change in the mass 
of the particle. An exact expression for the fidelity, appropriate for this perturbation, leads to an 
analytical result valid over very long time intervals, inversely proportional to the size of the mass 
perturbation. For such perturbations, we then calculate the long-time echo for semi-classical wave 
packets on periodic orbits. This paper also corrects an earlier calculation for a quantum echo, 
included in a previous paper, (nlin. CD/0307025 I. We explain the reasons for this correction. 



I. INTRODUCTION 

The search for signatures of classical chaos in corre- 
sponding quantum systems is one of the main themes of 
quantum chaos. The literature in this subject is large, 
much of it described in recent books by Haake and by 
Stockmann 0, Q • Here we consider a simple version of 
this question, namely, we describe a semi-classical calcu- 
lation of the short time spreading of a wave-packet for 
the quantum version of a classically chaotic system and 
show, as one might expect, that it is governed to a large 
extent by the Lyapunov exponents characterizing the ex- 
ponential separation of close trajectories of the classical 
system. We consider a Lorentz gas, which consists of a 
particle, or a collection of non-interacting particles, mov- 
ing in d-dimensions among a collection of fixed scatterers, 
generally taken to be d-dimensional hard spheres, with fi- 
nite horizon. The case where the scatterers are centered 
at the vertices of a regular lattice is the Sinai billiard. 
Related work on this problem has been done by Wirzba 
3] , which will be mentioned in Section II. 

We consider the propagator for a semi-classical parti- 
cle moving among the scatterers. The de Broglie wave 
length of the moving particle is taken to be small com- 
pared both to the size of a scatterer and to the average 
distance between scatterers. The propagator is evaluated 
by semi-classical methods for time intervals where a num- 
ber of collisions take place. As long as the wave-packet 
remains small, its spreading with time is governed by a 
set of equations that are the quantum counterparts to the 
curvature equations of Sinai and co-workers jj] that de- 
termine the Lyapunov exponents and Kolmogorov-Sinai 
(KS) entropy for the classical system. We can then eas- 
ily relate the spreading of the wave-packet to the classical 
Lyapunov exponents. 



Next we apply this result to a calculation of the wave- 
packet auto-correlation function (the return probability) , 
C(t), defined by 



C{t) 



-itH/h 



10} 



(1) 



where H is the Hamiltonian of the system and |0) is an 
initial quantum state. We find that, for wave packets on 
periodic orbits, this function exhibits a series of maxima, 
with amplitudes decreasing mainly exponentially with 
time, as exp(— /iks^)j where Hks — J2i ^ i s the classi- 
cal KS-entropy, and the Xi are the positive Lyapunov 
exponents for the corresponding periodic orbits. This 
type of decay of the auto-correlation function was first 
described by Heller [5|. Here we also calculate the coef- 
ficient in front of the exponential. This coefficient has a 
sub-exponential time dependence. The auto-correlation 
function exhibits a set of maxima separated by deep min- 
ima that appear for simple physical reasons, as we explain 
in Section III. 

These results are limited in applicability to times less 
than the Ehrenfest time, which is the time necessary for a 
wave packet to spread to the size of a scatterer. However, 
there is an application of them to the Loschmidt echo, 
or quantum fidelity, of a special type which is valid for 
a much larger time interval, greater than the Ehrenfest 
time. The Loschmidt echo 0,1113, M(t), is defined by 



M(t) = (0\e u ^ +s y h e- itH ^\0) 



(2) 



Here H is the Hamiltonian for the system, E is a small 
perturbing Hamiltonian, and |0) is some initial quantum 
state. For the case that H is the Lorentz gas Hamil- 
tonian with hard sphere scatterers, and the perturba- 
tion is a small change in the mass of the moving parti- 
cle, it is straightforward to show that M{t) is equal to 
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the wave-packet autocorrelation function evaluated at a 
scaled time, which can be made to be much shorter than 
the physical time t, by choosing a suitably small mass 
perturbation. Therefore for this special perturbation and 
hard sphere Lorentz gas system, the quantum fidelity can 
be evaluated for very long times, if one knows the behav- 
ior of the auto-correlation function for a much shorter 
time interval. 

An earlier version of this paper [To| contained the claim 
that the Loschmidt echo, defined above, decays over some 
time interval as e~ 2At . This claim is wrong, and our 
reasons for reaching this conclusion will be discussed in 
Section III. 

This paper is organized as follows: In Section II, we 
construct the semi-classical propagator for the moving 
particle and show that when it is used to determine the 
time evolution of a small initial Gaussian wave packet, 
the spreading of the wave packet is, for times less than 
the Ehrenfest time, determined by the classical curva- 
ture equations 0. These equations describe the rate of 
spreading of a classical, infinitesimal trajectory bundle. 
In Section III, we apply this result to calculate the wave- 
packet auto-correlation function for periodic orbits, as 
an illustration of the behavior predicted and described by 
Heller . In Section IV we derive an exact expression for 
the Loschmidt echo, M(t), for a quantum particle mov- 
ing in a hard sphere Lorentz gas, where the perturbed 
Hamiltonian is obtained from the unperturbed one by a 
small change of the mass of the moving particle. This 
result allows us to describe the behavior of the echo for 
long physical times in terms of the short time behavior of 
the wave packet auto-correlation function, C(t). We also 
discuss the error in our previous paper LA. The calcu- 
lations in these sections are for two-dimensional systems 
with hard disk scatterers. The three dimensional ver- 
sion of this work is presented in Section V. Here we show 
that the role of the positive Lyapunov exponent in our 
calculations for two dimensional systems is taken by KS 
entropy, i.e. by the sum of the two positive Lyapunov ex- 
ponents for the three dimensional system. We summarize 
our results in Section VI. 



II. THE SEMI-CLASSICAL PROPAGATOR FOR 
THE LORENTZ GAS 

We consider the semi-classical motion in two dimen- 
sions, d = 2, of a Gaussian wave packet, with average 
momentum p , whose initial form is given by 



(r|0) = ip (r) = (27rcr|| cr ) 



-1/2 



x exp [ — C 



T 



(3) 



where X — 7i/|po| is the de Broglie wavelength of the 
moving particle, a 2 Q = lfRe(il7*) and a 2 = l/Tle^ 1 ) 
characterize the size of the wave packet in the £- and re- 
directions respectively ( Tie denotes the real part). The 




FIG. 1: Particle-fixed frame of reference at time t — 0. 

(£, rf) system of coordinates is chosen with its origin at 
the center of the wave packet, r , and £-axis pointing in 
the direction of po, with ry-axis perpendicular to pp: 



r = r + U 



(4) 



where Uo is a 2 x 2 real matrix relating the two coordinate 
systems, see Fig. LI 

When the wave packet is far from any scatterers, its 
time propagation is dominated by free streaming, de- 
scribed by the propagator 



Gfs(r ' r ' t)= to) eXP 2^ (r - r) 



(5) 



where m is the mass of the moving particle, and d = 2. 
Application of this propagator to the wave function given 
by Eq. yields, up to an irrelevant phase factor, a new 
Gaussian wave packet of the form of Eq. J2J with 



f2|| — > Q» t = fi|i + iXvt/2 
fl Q -> n t = Oo + iXvt/2, 



(6) 
(7) 



where v = \po\/m is the average velocity of the particle. 
The new particle-fixed frame of reference is related to 
the stationary one by means of Eq. Q, with the wave 
packet center, r , replaced by r t = r + (po/\po\)vt and 
Ut = Uo- The average momentum of the wave packet 
stays unaffected: pt = Po- 

To find the semi-classical propagator describing a col- 
lision of the particle with one scatterer, we start with the 
general expression for the semi-classical propagator as a 
sum of terms of the form 



G S c(r,r',i) 



/ i \ d / 2 

I \ | D |l/2 



\2Ttih ) 



x exp i 



S(r,r',t) 7r/j 



(8) 



where S(r,r',t) is the classical action along a classical 
path from r' to r in time t, is an index equal to twice 
the number of collisions of the particle with hard disk 
scatterers over time t [12, D = det(-<9 2 Sy<9rcV), and 
d = 2. In general, there are two classical paths connect- 
ing points r and r', assuming that r is not in the geo- 
metric shadow of r': a reflected path and a direct one. 
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The contribution of the direct path from r' to r to the 
time evolution of the wave packet is negligible after time 
t if a classical particle with momentum po would collide 
with the scatterer during the interval (0,t). Thus, we 
only consider the propagator given by the reflected path. 
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FIG. 2: Particle-fixed frames of reference: (C',?/) a ^ time 
and (£, rj) at t. 

Consider a wave packet centered around r at time t = 
before a collision and around ro at t after the collision, 
Fig. |21 The system of reference originates at the point 
of classical collision. We suppose that the wave-packet 
to which the propagator will be applied is sufficiently 
small that we only need to find classical trajectories by 
minimizing the action for points starting close to rg and 
ending close to ro at time t. We then write the action, 
Sr^o + ^ r ) r o + <5r', i), for a trajectory from r + 5r' , 
\Sr'\ <C \r' \, colliding with a scatterer at point R and 
arriving at ro + Sr, \Sr\ <C |ro|, at time, t, as 



SrOo + Sr,r' 



It Vl 



5r',t) 



ro 



Sr - R 



(9) 



The variation of this action with respect to the point 
of collision, R, leads to an extremum equation that is 
used to determine the collision point, R. The algebra 
simplifies a bit if we take the case where \r' \ = \tq\ = r. 
Using the variational procedure, we find 

G SC (C, 77, t',rt', t) = G« (C, (', t)Gi 2 J (r), t/, t), (10) 



with 



/ 1 



\47riA"7 



1/2 7(C + 2r-Q 2 

exp — (11 
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and 



1/2 



x exp ■ 



4:iriXr(r + a cose; 
70(77 — r)') 2 cos (f> + 2ir{j] 2 + {rj') 2 ) 
AXr(r + a cose/)) 



(12) 



Here 77, £ and 77', £' are coordinates with origins at ro and 
r' Q respectively, such that (, £' are along the direction of 
the probability current, and 77, 77' are in directions per- 
pendicular to CjC'i respectively, as illustrated in Fig. El 
Further, <fi is the angle of incidence in the collision. The 
time dependence in the propagator appears in r, through 
the relation 2r — vt. 

There are limits to the range of applicability of the 
semi-classical propagator given by Eq. (fTU^ ■ First, the 
particle's wave function is supposed to be confined to a 
small region in space, the linear size of which is much 
smaller that the radius of the scatterer, throughout the 
time interval (0, t). It is this limitation that allows one to 
consider only the reflected path while deriving the prop- 
agator. Second, the wave packet size is assumed to be 
much smaller than the distance r between the center of 
the wave packet and the point where the particle would, 
classically, collide with the scatterer. This assumption 
makes it possible to expand the coordinates of points 
connected by the propagator about corresponding wave 
packet centers. 

The propagator in the direction of motion given by 
Eq. I|ll|) is simply the free streaming expressed in 
particle-fixed coordinate frames, showing that time evo- 
lution of ^-component of the time dependent wave packet, 
ipt, is unaffected by scattering events. The 77-component 
of the propagator, Eq. (|T2^I . can be easily shown to satisfy 
the identity 



(13) 



GH\ V , V ',t)= I d m I d V2 G ls (v,V2,t/2) 

xC(772,77i)G fs (77 1 ,77',t/2), 



where we introduced an instantaneous collision propaga- 
tor, (7(77,7/), according to 



C(t7j V ) — fi(v ~ v) ex P 



X a cos i 



(14) 



Eq. (|13fl allows to represent the propagator for a sin- 
gle scattering event, G sc (£, 77, r]', t), as a product of 
three successive propagators: (i) a free streaming prop- 
agator, Gf s (Ci, 771, r]',t/2), (ii) an instantaneous col- 
lision propagator, 5(^2 — Ci)G(?72, 77i) affecting the 77- 
component of ipt , and (iii) another free streaming propa- 
gator, G fs (C,77,C2,772, t/2). 

Assuming that the wave packet size at remains smaller 
than radius a of a scatterer, over the time t, we now con- 
struct the propagator for a trajectory with several col- 
lisions of the moving particle with scatterers as a com- 
bination of free particle and single collision propagators. 
This is appropriate in the semi-classical approximation 
when the size of the wave packet is small compared to the 
size of a scatterer, and to the average separation of the 
scatterers. Both free flight and instantaneous collision 
propagators leave the Gaussian form of a wave packet 
invariant. While the effect of the free streaming is de- 
scribed by Eqs. JHIIZIlj the instantaneous collision propa- 
gator, Eq. I|14f) . when applied to a Gaussian wave packet 
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leads to an instantaneous change in Q given by 



4/ 



fi- 



n- 



Xa cos i 



(15) 



where superscripts ± are used to distinguish variables 
immediately before and immediately after a collision. As 
mentioned above, is unaffected by instantaneous col- 
lisions: f2jj~ = Ojj~. 

The free streaming transformation of f2 t , coupled with 
the collisional transformation of f2~ to f2 + given above 
provides a direct connection between this semi-classical 
analysis of wave packet motion and the method of Sinai 
et al. for analyzing the ergodic properties of the classi- 
cal Lorentz gas in terms of the curvature of a classical 
wave front 0, Il3| . In fact a simple transformation al- 
lows us to recover the classical equations, and to identify 
the appearance of the positive Lyapunov exponent in the 
semi-classical formulae. To see this let us define complex 
radii of curvature, p\\ and p, according to 



fin = — Xpn and = — Xp. 

ii 2 n\ 2 v 



In terms of p, Eqs. iffl and H~5|l read 



Pt = Po + vt 
1 1 2 



free streaming, 
collision, 



(16) 

(17) 
(18) 



while p\\ t = p|i + vt regardless of whether any scattering 
events have taken place over time t. These equations 
for p are identical with the curvature equations 4] for 
the classical Lorentz gas. In an unpublished manuscript 
describing the diffractive scattering of a wave packet by 
a circular scatterer, Wirzba Q noted that the curvature 
equations can also be extracted from his formalism. 

To describe the spreading of a Gaussian wave packet 
in the Lorentz gas, we consider a sequence of collisions 
parameterized by a set of times {tj} together with a 
set of collision angles {cf>j}. Direct substitution of the 
free streaming transformation for p t , Eq. 117|) . into the 
expression for the size of the wave packet along the r\- 
coordinate, i.e. along the direction perpendicular to the 
average momentum p of the particle, of = 1/ Tie^l^ 1 ) = 
X/(2Jm(pi 1 )), yields 



<?t = <?t, 



hi +v{t-t j ) 



Of, exp v Tie 




for U <t <t 



j+l- 



(19) 

It follows from the relation between a 
and p, and the change in p on collision, that the instanta- 
neous scattering transformation does not change the size 
of the wave packet (af. = o t ~). Thus, we can propagate 
at backward in time to get 



a t = <To exp [ v lie \ — 

o Pt' 



„*A t 



(20) 



where oo is the initial size of the wave packet at t = 0, 
and 



At = -r 7& / — = - / — 



Pr t J p 



(21) 



where we introduce a real radius of curvature, p, as 

l/p=fie(l/#. (22) 

The quantity At can be thought of as a wave packet 
stretching exponent over a time t. It differs from the 
classical Lyapunov exponent A because it contains quan- 
tum effects and the limit of infinite time is not taken. 
The stretching exponent, At, converges to the Lyapunov 
exponent, A, in the long time classical limit: 



lim lim At = A. 



(23) 



In order to prove Eq. I|23|l . one needs to show that p 
becomes the classical radius of curvature for the classi- 
cal Lorentz gas as A — > 0. Substituting Eq. 12211 along 
with Xm(\jp) — X / {2a 2 ) into the transformations for p t , 
Eqs. one gets 



Pi 



(go + vt) 2 + e {vtf 
Po + vt + e vt 



a t = —^{p + vt) 2 + e {vt) 2 
Po 



free streaming, 



(24) 



and l/p + = l/p + 2/(acos</>) together with a + = a 
at a collision. Here 



Xp t 
2a 2 



(25) 



contains all the quantum effects; it vanishes as X — > 0, 
which makes Eq. 124|) converge to its classical counter- 
part |J, Il3| . Another way to visualize the semi-classical 
corrections is to rewrite Eq. (|24|) in differential form: 



pt = v(l - Et) and a t = va t /p t - 



(26) 



Here the second equation has its classical form, and the 
quantum correction is apparent in the first equation: it 
shows that the free flight spreading of the wave packet 
results from a combination of a classical linear separation 
of trajectories and the quantum spreading due to the 
Uncertainty Principle. 

The role of the Uncertainty Principle becomes appar- 
ent from the following simple consideration. Suppose one 
prepares a tiny minimal wave packet with spatial uncer- 
tainty <tq. The corresponding uncertainty in momentum, 
Ap, is then given by ooAp = h/2. After some time t the 
wave packet size evolves to ctup ~ (Ap/m)t — Xvt/(2ao) 
merely due to the Uncertainty Principle. Writing the ge- 
ometrical (classical) spreading as ctcl = o-q(1 + vt/ po), 
we notice that at in Eq. (|24J) is essentially a simple com- 
bination of a ql and oup , namely Of = \fo\ 



'cl 



'up- 
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FIG. 3: Free flight time evolution of p and er: (a) classical 
case, e = 0, (b) quantum case, e > 0. 



Fig. |3 illustrates the free flight dynamics of p and a 
given by Eq. 1(21)1) . Fig. |3fa) pictures the classical limit, 
A = 0: an arc of instantaneous radius p and length a 
moves with constant velocity v along the "cone" orig- 
inating at a point O. Eq. (|26|l with e t = describes 
the time evolution of p and a in this case. In quantum 
regime, e > 0, the point O is also moving in the same 
direction as the arc, but with a different, time depen- 
dent, velocity equal to ev, see Fig.[3Jb). It can be shown 
from Eqs. I|25l I26j) that e t ~ t~ 2 as t — > oo, implying the 
convergence of point O to some point O' in the long time 
limit, see Fig.[3fb). The time evolution of p and a is then 
dominated by classical equations when O is close to O'. 

We now define an interval of time, called the Lyapunov 
regime for which the values of p and a satisfy the inequal- 
ity e <C 1. It follows from the free flight and collision 
transformations for p t and a t that et is a rapidly decreas- 
ing function of time, see Fig. 0] Therefore, once in the 
Lyapunov regime the system stays in it for some time ii, 
at which a becomes comparable with the size of scatter- 
ers, and our collision analysis breaks down. It can also be 
shown that if the Lyapunov regime inequality is not sat- 
isfied at t — 0, and the wave packet is small, the system 
rapidly evolves to a state for which this inequality is sat- 
isfied. During this transient regime p rapidly decreases 
whereas a does not change significantly, see Fig. 0] 

In the Lyapunov regime, the first equation in Eqs. I|26|) 
reduces to its classical counterpart, p « v, so that a t 
grows exponentially in the same way as a small pencil of 
trajectories separates exponentially in the classical sys- 
tem. That is at — <Jq exp (t\t), where At is given by 
Eq. H21|) and calculated using only classical mechanics. 
It is useful to remark that At typically reaches a value 
close to the classical Lyapunov exponent A after only a 
few collisions, see Fig.0] On the other hand, one can es- 
timate that the upper limit of th, the maximum duration 
of the Lyapunov regime, is t™ ax ~ (1/(2A)) ln(a/A), that 
is, about half of the Ehrenfest time which for suffi- 
ciently small A can be long enough for the wave packet 
to exhibit exponential spreading. 

Finally, we illustrate the exponential spreading of a 
Gaussian wave packet for the case of particles moving in 
short, periodic orbits. We numerically evaluate a t and pt 



in Eq. Il26[l (and in the p~ to p + collision transformation) 
for the simplest periodic orbit: a particle moving back 
and forth along the line connecting the centers of two 
disks. Fig. 21 shows at, pt, and quantity et, given by 
Eq. i|25|) , for the two disks of radius a = 1 , and the center- 
to-center separation R = 3. The particle is placed in the 
middle between the two disks at t = 0, and has the de 
Broglie wavelength A = 10~ 7 . The initial wave packet 
is characterized by uq = 2 • 10~ 4 and po = 10, so that 
£o ~ 156 and the system is far from the Lyapunov regime 
at t = 0. Fig.0|shows that it only takes a single collision 
for the system to reach the Lyapunov regime, e <C 1. 




FIG. 4: Wave packet size, a, real radius of curvature, p, and 
e, are shown as functions of time, t, for a two-disk periodic 
orbit. Disk radii a = 1, center-to-center separation R — 3, de 
Broglie wave length X = 1CP 7 . The corresponding Lyapunov 
exponent X/v ~ 1.32. Initial wave packet size ao = 2 • 10 -4 
and po ~ 10. The particle is located in the middle between 
the two disks at t = 0. Exponential trends are shown for plots 
of a and e. All distances are measured in units of a. 

The parameters in Fig. 0] are chosen so as to illustrate 
the essential regimes: a short decay of quantum effects 
(e becomes less than unity), followed by the Lyapunov 
spreading of the wave packet, a ~ exp(Ai). The rela- 
tively small value of the de Broglie wavelength used in 
this example can be indeed achieved experimentally |14| . 

The classical Lyapunov exponent of a two-disk periodic 
orbit is known exactly Eg > 



R-2a 



ln R^±vmEM. (27) 



In our case Eq. H27J1 gives X/v w 1.32. The numerical 
evaluation presented in Fig.0]shows that a single collision 
is enough to initiate the exponential growth of the wave 
packet (with the rate given by the classical Lyapunov 
exponent), which persists for about 5-6 collisions. Our 
results do not apply for times longer than the duration 
of the Lyapunov regime, vt^/a « 6. 
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III. THE WAVE-PACKET 
AUTO-CORRELATION FUNCTION 



As an application of the analysis of wave packet dy- 
namics developed above, we calculate the wave packet 
auto-correlation function, C(t), defined in Eq. (JTJ, for 
particles moving in periodic orbits. Here the initial state, 
|0), describes a Gaussian wave packet centered about r 
with its average momentum p , such that the phase point 
(r , po) lies on a periodic orbit of the corresponding clas- 
sical system. The auto-correlation function for periodic 
orbits in billiard systems was studied by Heller y| some 
time ago using different techniques. The calculations pre- 
sented here agree with Heller's results and provide some 
additional information about this correlation function. 

The reasons for restricting our calculations to periodic 
orbits are as follows. The expansion used above to obtain 
the semi-classical single collision propagator in the pre- 
vious section, Eq. ijlfjl) . is correct for wave packets which 
are small compared to disk radii and average separation 
among scatterers. Mathematically, this limitation is a 
consequence of the truncation of the expansion of the 
coordinates of starting and final points connected by the 
propagator, r' and r respectively, about the centers of ini- 
tial and final wave packets respectively, i.e. r' = r + Sr' 
and r = ro + Sr. Therefore, one gets a close approxima- 
tion to the particle's wave function at positions close to 
the wave packet center, ro, but the approximation may 
fail on the periphery of the wave packet. Our calculations 
of the auto-correlation function, C(t), or Loschmidt echo, 
M(t) 7 are only reliable when the relevant overlap inte- 
grals are dominated by central region of the wave func- 
tion, and contributions coming from wave packet wings 
can be neglected. This condition is most easily satisfied 
when the classical motion is along a periodic orbit. The 
result for the decay of the Loschmidt echo reported earlier 
[Tot . as exp(-2Ai), is incorrect because this restriction 
on the validity of the semi-classical approximation was 
not properly taken into account. This error first became 
apparent when we obtained an exact result that is incon- 
sistent with this decay. The exact result is described in 
the next section. 

Consider a wave packet whose initial average coordi- 
nate, ro, and momentum, po, correspond to a phase space 
point on a periodic orbit, with period T, of the classical 
Lorentz gas. Suppose that nT is smaller than the Ehren- 
fest time, for n = 0, 1,2, ...,N so that we can apply 
the analysis developed in the previous section in order to 
propagate the wave packet over times 



t = nT + A/v, 



(28) 



where the displacement A is sufficiently small in order for 
the initial and final wave packets to overlap significantly, 
as illustrated in Fig.[S] For simplicity we take the initial 
wave function, ipo, to be a "circular" wave packet, i.e. 




FIG. 5: A Gaussian wave packet is shown at time t = 
(centered about point O), and at a later time t = nT + A/v 
(centered about point P). Points O and P lie on the same pe- 
riodic orbit, and are separated in time by n (= 0, 1, 2, . . . , iV) 
periods, T, of the periodic orbit, plus a short time interval 
A/v. The separation distance A is assumed to be sufficiently 
small in order for the initial and final wave packets to overlap 
significantly. 



— and pyo = pa with the explicit form, 



x exp 



/ 1 



V27TO-2 



1/2 







2i \ 






*Po) 



(x 2 +y 2 ) + -x 



(29) 



where x-axis is directed along po. In the same coordinate 
system the wave packet ^ t propagated from tpo over time 
t, given by Eq. H28|) . reads, up to an irrelevant phase 
factor, 



1 



27TCT|| cr 

-(- 

4 U 2 



1/2 



exp 



4 W *pii/ v 



2i 
~X~p 



(30) 



Here, an , a, pn and p depend on time t through a se- 
quence of free flight and collision transformations devel- 
oped in the previous section. The probability distribu- 
tion | V^o C r ) 1 2 i s negligible outside a small circle of radius 
r ~ do- Therefore, the main contribution to the over- 
lap (ipolifit) comes from the points inside this circle, and 
the central regions of the wave packets will dominate the 
integrals for small center-to-center separations A, as il- 
lustrated in Fig. O 

A straightforward calculation shows that for t given by 
Eq. J2HI) 



C(t) 



dr ip*(r)ip t (r) 



-exp(-cvA 2 ), (31) 
a 
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where 
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A = -x 



-Re 



9\\9 

9\\ 




(32) 
(33) 



with 





' 1 


9 = 






J* 




1 


9\\ - 







2i 



1 1 

P Po 
2i ( 1 1 



(34) 



As seen from Eq. (|31|l the auto-correlation function ex- 
hibits a sequence of peaks corresponding to partial re- 
construction of the wave packet at times t — nT. These 
peaks, first studied by Heller have a simple phys- 
ical origin: the wave packet repeatedly passes through 
the starting point giving rise to strong maxima of the 
return probability C(t). These maxima are periodic or- 
bit revivals and should be distinguished from more gen- 
eral classes of quantum revivals that do not require a 
particular periodic orbit for their appearance [||. It 
can be shown that time dependence of A and a is sub- 
exponential compared to to the exponential growth of a 
with time, so that the periodic orbit revival peaks have 
predominantly Gaussian form. The strength of the peaks 
decreases exponentially with time with a rate given by 
the Lyapunov exponent of the periodic orbit. This fol- 
lows from the fact that height of the peaks is mainly 
determined by the exponential growth of the size of a 
wave packet, a ~ exp(Ai). It is worth noting that A 
decays with time in a power law manner making the 
auto-correlation function to decay slightly faster than 
exp(— Xt), see Fig. El 

Fig. shows the numerical evaluation of the revivals 
in Eq. I|31|l for the two-disk periodic orbit described in 
previous section, see Fig. A particle of the de Broglie 
wavelength X = 1CP 7 moves back and forth between two 
disks of radii a = 1, with the center-to-center separa- 
tion R = 3, along the line connecting the centers. The 
initial wave packet is located in the middle between the 
two disks, and is characterized by cr|| = Co = 2 • 1CP 4 
and p|| = po = 10. The left part of Fig. shows the 
revival maxima C max , which occur at i ma x = nT. The 
right part shows the auto-correlation function in small 
neighborhoods of the corresponding maxima. 

It easy to show that in case of the two-disk periodic or- 
bit, the revival peaks are separated by deep minima of the 
auto-correlation function. The minima occur when the 
average momenta of the original wave packet and the one 
propagated in time t are pointing in opposite directions. 
When this happens, ipo and ipt interfere destructively, 
and the auto-correlation integral, C(t), is very small. We 
calculate the overlap of -00 an d ipt at t = (n + 1/2)T, 
when the two wave packets are centered about the same 
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FIG. 6: Revival peaks of the wave packet auto-correlation 
function, C(t), for the same two-disk periodic orbit as in 
Fig. |1| a = 1 and R = 3. Particles de Broglie wave 
length X = 10 -7 . The initial wave packet is characterized 
by (To = (Tuo = 2 • 10~ 4 and po = p«o = 10. The exponential 
trend is indicated by a straight line. 



point, but move in opposite directions. The initial wave 
function is given by Eq. (gUJ, while 



Then, 



1 



27r<7||(7 



1/2 



exp 

1 

4 




C{t) 



A 



■ exp 



X 



2 T^9\\ 



(36) 

where A and gt\ are defined in Eqs. (|32l I34f> . It can be 
shown that the exponential in Eq. (|36|l is a very small 
number if the condition X <C (Tq,Pq is satisfied, e.g. in 
case of the periodic orbit considered above (Fig. Q] and 
Fig. 01 Veg\\ ranges from 2 • 10~ 8 to 4- 10~ 8 making C(t) 
smaller than (A/a) exp(— 4-10 6 ), which is practically zero 
when compared with values of C(t) at the PO revival 
maxima. 

Similar deep minima of the auto-correlation function 
will occur for particles in more complicated periodic or- 
bits, making the revivals very pronounced. 



IV. A SIMPLE HARD SPHERE LOSCHMIDT 
ECHO 

Having described the wave-packet auto-correlation 
function for the Lorentz gas, we now present a simple 
identity that allows us to use this correlation function 
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to calculate, analytically, the Loschmidt echo for a very 
particular perturbation. The Loschmidt echo was defined 
by Eq. We suppose that the perturbed Hamiltonian 
is obtained from the unperturbed one by changing the 
mass of the moving particle, from to to to + 5m. The 
identity depends on the fact that for hard scatterers, no 
matter what their shape, the eigenfunctions depend on 
wave numbers rather than on the mass of the moving par- 
ticle. That is, the sets of eigenfunctions for particles of 
different masses are the same, only the values of the en- 
ergy corresponding to the same wave numbers differ. The 
wave functions are the solutions of the scalar Hclmholtz 
equation 



(V 2 + k 2 ) k (r) = 



(37) 



which satisfy the Dirichlet boundary condition that 0k (r) 
vanishes on the surface of each scatterer, and on the 
boundaries of the system. 

We can express the time propagator for a moving par- 
ticle of mass to under a Hamiltonian H m as 



exp 



it 
n 



<)exp ( -^-k 2 
2m 



(0k |, (38) 



where the summation is over all possible eigenstates of 
the system. These eigenstates, |0k), satisfy an orthonor- 
mality relation 



(0k' 1 0k) = <5k' 



k',k, 



(39) 



times can be expressed in terms of the short time auto- 
correlation function. 

The physical origin of this result is straightforward. 
Classically, the perturbed and unperturbed masses follow 
exactly the same trajectory, but with different velocities. 
Hence the forward motion with mass m followed by the 
reversed motion with mass m + 5m has a final position 
which is different from the initial position, and corre- 
sponds to motion over the part of the path that is not 
reached by the time reversed trajectory. This is exactly 
reflected by the operator identity, Eq. I|40l) . Perhaps the 
most remarkable thing about this result is the fact that 
although for long times a small wave packet will spread 
over large distances, the Loschmidt echo in this case is de- 
termined by the short time spreading of the wave packet, 
even if the physical time is quite large. As mentioned 
in Section III, the exact result, Eq. I|42|l . is inconsistent 
with the Lyapunov decay reported earlier [Ioj|. 

As discussed in Section III, the auto-correlation func- 
tion C(t) exhibits a sequence of sharp revival maxima 
when the particle moves on a classically periodic orbit. 
The maxima occur at times i max multiple of the period T 
of the periodic orbit, and C(i max ) ~ exp(— At max ), where 
A is the corresponding classical Lyapunov exponent. Ac- 
cording to Eq. (|42[1 the mass perturbation Loschmidt 
echo M at time t is simply the auto-correlation func- 
tion C at the scaled time t s given by Eq. (|41|l . Thus like 
the auto-correlation function, the Loschmidt echo M{t) 
exhibits a periodic sequence of maxima at times 



where the choice between Kronecker and Dirac delta 
functions is dictated by the nature of the eigenstates. 
Eqs. (|38l I39fl hold for systems with hard wall potentials 
in any number of spatial dimensions. 

This representation of the time displacement opera- 
tor, Eq. I|38(l . together with the orthonormality condition, 
Eq. I|39|) . leads to the following identity 



t' = 

''max 



TO + 5m 

5m 



-U 



to + 5m 
n — -z T, 

dm 



(43) 



where n — 1, 2, . . . , N, such that NT is smaller than the 
duration of the Lyapunov regime, t^. The envelope of 
the maxima exhibits a mainly exponential decay: 



< M> 1 jrH m+Sm j <•:<]) 



it 



H m = exp 



it s 



-H„ 



where t s is a scaled time, related to the physical time t 
by 



t. 



5m 



-t. 



(41) 



M (Cax) = C(*max) ~ exp(-At max ) = CXp(-A s t max ). 

(44) 

Here we introduced a scaled Lyapunov exponent accord- 
ing to 



5m 



to + 5m 



-A. 



(45) 



This identity permits us to express the Loschmidt echo, 
for this special perturbation, in terms of the wave-packet 
auto-correlation function as 



M(t) 



i "it 

(0| exp ( —H m+ s 



(0| exp 



exp 



it 
n 



10} 



it. 



H„ 



10) 



(42) 



C(t s ). 



This is the main result of this section. For small per- 
turbations, 5m/ m <C 1, the Loschmidt echo for long 



It is important to note that the behavior of the Loschmidt 
echo described by Eq. (|44|) can persist for times much 
longer than t^ (for sufficiently small 5m) despite the fact 
that the analysis of the wave packet dynamics presented 
in Section II is valid only for times shorter than t^- 

The Hamiltonian perturbation used in this section is 
rather trivial since the perturbed Hamiltonian commutes 
with the unperturbed Hamiltonian. Therefore the results 
of this section are not to be compared with those obtained 
for more complicated perturbations such as distortion of 
the mass tensor 0, || . 
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GENERALIZATION TO THREE 
DIMENSIONS 



(& > ^1 j m ) an d (C2 j % j V-P ) J ust before and after the 
collision respectively, reads 



The derivation of the wave packet propagator in Sec- 
tion II, and the calculation of the periodic orbit revivals, 
Section III, were carried out for hard-disk systems in two 
dimensions, d = 2. We now generalize these calculations 
to the three dimensional case, d = 3, using methods sim- 
ilar to those used to describe the classical separation of 
close trajectories 0, 0| 

The initial Gaussian wave packet in three dimensions 
reads 



(r|0) = Mr) = 



3/4 



1 



a \\0 a 



4n,| 



1/2 



(46) 



where C^-axis is directed along the momentum p , see 
Fig. ^ while q = (t/ 1 ) , ?/ 2 )) T lies in the plane perpen- 
dicular to po; fio is a 2 x 2 complex symmetric ma- 
trix, and T-superscript denotes transposition. As in 
two-dimensional case, the origin of the orthogonal frame 
(CiV^iV ) travels with the center of the wavepacket 
with fixed axes, except at collisions, when the axes ro- 
tate so that the new £ axis is in the direction of motion 
of the center of the wave packet, see Fig. ^ The wave 
packet size in ^-direction a^ Q = 1/ 7?e (rijjg 1 ), while in 
rj-plane 



det Hein^ 1 ) 



Application of the free streaming 
Gf s (r,r',i), given by Eq. J5J with d = 
wave function above changes J7o to 

n t = n + i(Xvt/2)i, 



(47) 



propagator 
3, to the 



(48) 



where 1 is the 2x2 unit matrix; the change of the £- 
directional component of the wave packet is the same as 
in the two-dimensional case, Eq. (|SJ. 

The single-sphere scattering propagator is given by 
Eq. © with d = 3. As in the two-dimensional problem, 
only the reflected path contributes to the propagator for 
a wave packet small compared to the sphere radius, a. 
Closely following the arguments of Section II in three di- 
mensions, one can verify that the scattering propagator 
G sc (r, r', t) can be written as 



G sc (r,r',i) = J dr x J dr 2 G fs (r, r 2 , t/2) 



(49) 



xC(r 2! r 1 )G fe (r 1 ,r / ,t/2), 



where, in order to simplify the algebra, we consider the 
case that the corresponding classical collision takes place 
at time t/2. The instantaneous collision transformation 
G, when expressed in particle- fixed coordinate frames 



1 M (5°) 



x exp — r il Q(0,6»)?7 1 , 
Xa 



where 



and 



Q(4>,9) =P e diag 



- , COS ( 



cos 9 
sin 6 



- sin 9 
cos 9 



(51) 



(52) 



Here <fi is the angle of incidence in the collision plane, see 
Fig. 21 and 9 is the azimuthal angle that j^^-axis makes 
with the collision plane. Note, that the coordinate frames 

(Ci: 7 7i 1 \ 7 7i 2 ' ) ) an d (C2, W^^vi ) are related to each other 
by the 3x3 reflection matrix (I3 — 2nn), where I3 is the 
3x3 unit matrix, and n stands for the three-dimensional 
collision vector, as illustrated in Fig. [21 

As seen from Eq. (|5U|) the instantaneous collision does 
not affect but changes the ^-component of the wave 
packet according to 



n -K+) = n-K-) - il Q ( 

Xa 



(53) 



Introducing the radius of curvature matrix p as 11 = 
iXp/2 we obtain the three-dimensional equivalent of 

Eqs. (jninsji, 



p t = p + vtl free streaming, (54) 
) collision. (55) 



9 

-Q( 

a 



Both transformations preserve the symmetry of the com- 
plex matrix p. 

As in two-dimensional case, we consider a sequence of 
collisions parameterized by a set of times {tj} together 
with a set of collision angles {(j)j,9j}. Substitution of 
the free streaming transformation, Eq. Ij54(l . into the ex- 
pression for the size of the wave packet in the r/-plane, 



a\ r = X/(2J det Tmp t ), yields 



det 



Pu + v ( t - t j) 1 



det p t 



of. exp Tie J dr trp T . 



(56) 



for tj < t < tj + i. Here we used the identity det Tmp^ 1 = 
I det p\~ 2 det Imp. By propagating at backward in time 
we find 

a\ = 0% exp {v TieJ^ dr trp^ 1 ^ = erg e tht , (57) 
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where trg characterizes the wave packet at t = 0, and 



VI. SUMMARY 



h t = - n 



/ dr txp T 
Jo 



(58) 



The quantity ht is a time-dependent stretching exponent, 
which describes growth of the area of wave packet cross 
section perpendicular to the direction of particle's mo- 
tion. It can be shown to converge in the long time classi- 
cal limit to the classical Kolmogorov-Sinai (KS) entropy 
^■KS) equal to the sum of all positive Lyapunov exponents 
in the system: 



lim lim ht = /iks = / <\p 

\j>0 



(59) 



To complete the analogy with the two-dimensional 
problem we define a real radius of curvature matrix p 
and a real 2x2 matrix S in accordance with 



y(SS T )~ 1 



(60) 



It is easy to show that S determines the size a of the wave 
packet, a 2 — | det S|, and is not affected by the collision 
transformation given by Eq. I|55|) . while p satisfies 



p- 1 ^ =p- 1 (-> + -Q( 



(61) 



at collisions. The free streaming time evolution of p and 
S is given by the differential equations 



1 dp 

— r = 1 

v dt 



T\-2 



1 rfS 

v dt 



(62) 



which are the three dimensional version of Eqs. (|26J) . 
Since S + = E~, the second equation in Eqs. H62J) can 
be integrated to get 



S t = T exp 



1 / drp T 
Jo 



So, 



(63) 



where T stands for the time ordering operator. Finally, 
taking the determinant of both sides of Eq. I|63|l we re- 
cover Eq. (|57|l . namely of = a 2 exp(th t ) with 



ht = - 



- / drtrp T . 
1 Jo 



(64) 



In this paper we have considered the short time spread- 
ing of a small Gaussian wave packet for a particle moving 
in an array of fixed, hard sphere scatterers, in both two 
and three dimensions. Our calculations are based upon 
the semi-classical expression for the quantum propagator 
in terms of the classical action for paths of the particle. 
We find that for times less than the Ehrenfest time, the 
spreading of the quantum wave packet is determined by 
the sum of the positive Lyapunov exponents that describe 
the classical separation of nearby trajectories. We used 
the expressions for the propagator to calculate the wave 
packet auto-correlation function for periodic orbits. Our 
results agree with earlier results of Heller [5j: (1) this 
function exhibits a set of sharp maxima, the periodic or- 
bit revivals, whenever the moving wave packet overlaps 
with the initial one and has the same velocity direction; 
and (2) The strengths of the maxima decrease exponen- 
tially with a decay rate given by the positive Lyapunov 
exponents. When the velocities are oppositely directed, 
the correlation function takes on extremely small values, 
even though the wave packets spatially overlap. Finally 
we used a special property of the eigenfunctions for hard 
sphere Lorentz gases to evaluate the quantum fidelity, 
or Loschmidt echo, for a perturbing Hamiltonian that is 
just a small change in the mass of the moving particle. 
The property that the eigenfunctions are independent of 
the mass of the particle, when expressed in terms of the 
wave number, allowed us to relate the Loschmidt echo at 
long times to the wave packet auto-correlation function 
at much shorter times. Therefore, for periodic orbits, 
at least, the Loschmidt echo will exhibit the same kind 
of periodic orbit revivals as one finds for the correlation 
functions. 

It would be very interesting if one could provide an- 
alytic calculations of quantum echos and revivals over 
longer time intervals for Lorentz gases or other, simpler 
models, such as quantum multibaker models [T3|. We 
would need to find appropriate techniques for analyzing 
the space and time development of wave packets for times 
longer than the Ehrenfest time. This would enable us to 
describe the numerical results for the Lyapunov decay 
obtained by Pastawski, Jalabert and co-workers d,EJ. It 
is not yet clear to us how this might be accomplished. 



We can also calculate the wave packet auto-correlation 
function, C(t), defined in Eq. Q, on periodic orbits for 
times t given by Eq. J5SJ. The algebra is straightfor- 
ward, but rather lengthy, and we provide here only the 
main result of the calculation: the auto-correlation func- 
tion, C(t), exhibits a sequence of sharp maxima, periodic 
orbit revivals, which occur at times t max — nT, with 



a t , so that the envelope of the PO revivals 



shows mainly exponential decay with the rate given by 
the KS-entropy, C max ~ exp(~h KS t 

max } ■ 
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